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Abstract. In lattice QCD, obtaining properties of heavy-light mesons has been easier 
said than done. Focusing on the B meson's decay constant, it is argued that towards the 
end of 1997 the last obstacles were removed, at least in the quenched approximation. 
These developments, which resulted from a fuller understanding and implementation 
of ideas in effective field theory, bode well for current studies of neutral meson mixing 
and of semileptonic decays. 
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Eleven years ago at Lattice '87, three talks gave birth to a new field of research, 
A ■ the study of heavy quarks in lattice QCD. Estia Eichten [1] emphasized that lat- 

tice QCD should provide reliable information on hadrons with heavy quarks, which 
Qh! would help determine the Cabibbo-Kobayashi-Maskawa (CKM) matrix. He sug- 

gested starting with the static approximation and adding corrections in I/tuq. One 
of the developments to grow from this suggestion is the heavy-quark effective the- 
ory (HQET), an enormous subject in its own right. Peter Lepage [2] introduced 
non-relativistic QCD (NRQCD) as a tool for studying spectroscopy and matrix 
elements of systems with one or more heavy quarks. This effective theory also 
has enjoyed widespread application. Finally, Luciano Maiani [3] presented, among 
other things, the first well-known calculation from lattice QCD of the leptonic de- 
cay constant of the B meson. Maiani and his collaborators took Wilson's action 
for light quarks and bravely applied it to the b quark, even though their lattice's 
ultraviolet cutoff was below the mass mj. 

The idea that lattice QCD could play a role in interpreting future experiments 
was exciting, and it attracted much attention. In addition to the spectra of 
quarkonium and "heavy-light" hadrons, considerable effort has been devoted to 
the hadronic matrix elements for leptonic and semileptonic decays of heavy-light 
mesons (B, B s , D, and D s ) and neutral meson mixing. The decay constant of the 
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generic heavy-light pseudoscalar meson, denoted fp, parameterizes the hadronic 
amplitude for leptonic decays. It has received the most attention, especially fp. 
It was expected to be the most straightforward of matrix elements and, thus, a 
bellwether. Unfortunately, a reliable calculation could not be done quickly. For- 
tunately, the technical and conceptual difficulties have been largely overcome, and 
the calculation of fp in lattice QCD has grown up at last. 

Computationally fp is indeed just as easy as f n , but the interpretation of the 
results has not been obvious. Consequently, the literature contains a wide range 
of estimates, some of which should not be taken seriously. For example, several 
early calculations in the static limit contained too much contamination from radial 
excitations B [4], and the results are misleadingly high. 

More recently, a fuller understanding of the interplay between the effective theo- 
ries and the lattice has helped to reduce the effects of lattice artifacts. The effective 
theories NRQCD and HQET are derived from QCD by lowering the renormaliza- 
tion point, or cutoff, /i until |p| <C fi ~ rag, where |p| is the heavy quark's typical 
momentum and tuq its mass. Because \p\/rriQ <C 1 the interactions in the effective 
Lagrangian can be organized in powers of \p\/rriQ. (In quarkonium \p\ ~ a s mq; in 
heavy-light systems |p| ~ Aqcd ~ 200 MeV.) The two effective theories share the 
same Lagrangian, although the power of \p\jmQ assigned to operators of higher 
dimension can differ. One can take the effective theories' Lagrangian and intro- 
duce the lattice as an ultraviolet regulator, choosing the lattice spacing a so that 
ttiq pd fi ~ ir /a. This lattice theory is often called lattice NRQCD [5]. There are 
discretization effects, which are the higher- dimension operators multiplied by cal- 
culable coefficients. Although some of these operators are new, many are the same 
as in the I/tuq expansion. Consequently, physical 1/ttiq effects and artificial a 
effects have become intertwined. Lattice practitioners must disentangle them and 
remove the lattice artifacts, at least to the desired accuracy. 

Alternatively, one can start with an action derived for m q a <C 1, such as Wilson's, 
and ask what happens when one applies it for rriQa ~ 1. This, essentially, is the 
way of Ref. [3]. Since then, many experts have said (and still do, out of habit), that 
a heavy quark cannot be put directly on the lattice because rriQa ~ 1. On the other 
hand, numerous calculations have been published with mqa ~ 1, so there must be 
more to the story. Indeed, the lattice theory does not break down. Instead, the 
lattice artifacts are again intertwined with the I/tuq effects, as in lattice NRQCD. 
The same operators appear, but the coefficients are different, though still calculable. 
For a class of actions based on the Wilson action, it has been shown, to all orders 
in perturbation theory, that the coefficients remain small, for all mqa [6]. 

In preparing a brief review of a subject one is faced with the choice between 
a catalog of all recent results or a synthesis of developments over a longer period 
of time. The proceedings of the Lattice conferences provide excellent examples 
of the former [7-9]. By contrast, this paper gives a view of the (theoretical and 
computational) progress, focusing on fp. Owing to space limitations the material 
presented on the allied subject of neutral meson mixing and on phenomenologically 
promising form factors of semileptonic decays is brief. 
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NUMERICAL LATTICE CALCULATIONS 



When one thinks of lattice QCD, one usually thinks of large-scale numerical 
calculations. This approach computes the functional integrals of quantum field 
theory by applying a Monte Carlo method with importance sampling. Statistical 
errors arise here, and with more and more computer time these errors can be made 
arbitrarily small. Over the years various clever techniques have been devised to 
enhance the "signal-to-noise" ratio, that is, to reduce the statistical error for fixed 
computing resources. 

To use Monte Carlo methods, three modifications are introduced. First, space- 
time becomes a finite box, usually with periodic boundary conditions. Second, the 
spacetime continuum becomes a discrete lattice. Last, the so-called quenched ap- 
proximation is applied. The first two are common to many kinds of numerical anal- 
ysis, but the last requires a short explanation. The quenched approximation treats 
a hadron's valence quarks and all the exchanged gluons fully, including retardation, 
but the back-reaction of closed quark loops on the gluons is omitted. For particle 
physics a more descriptive name (and one of the original names) would be the va- 
lence approximation, but the term "quenched," taken from statistical mechanics, is 
more commonly used. The back-reaction of closed quark loops is computationally 
burdensome. Because its omission saves computer time, the quenched approxima- 
tion is a useful way to control the other errors and, thus, to teach theorists how to 
analyze uncertainties. 

Let us return to the first two approximations. A volume larger than a few fm 
on a side should be good enough. After all, one does not expect the true size and 
boundary of the universe to effect the physics of hadrons. It is with the lattice 
itself that the subject becomes a craft. The continuum limit can be reached by 
taking the lattice spacing a — > with brute force, or by improving the action to 
reduce discretization effects, or by a combination of the two. The crudest form 
of brute force, namely to take rriQa <C 1, would require, for the bottom quark, a 
fantastically small lattice spacing. This has never been done. 



Table 1 contains a list of specific reactions with conventional parameterizations 
of hadronic matrix elements. Together each pair can determine the listed element 
of the CKM matrix. Every element of the CKM matrix appears except Vtb and the 
entire (3x3 unitary) matrix can be constrained. 

In the Standard Model, the leptonic partial width of a generic pseudoscalar 
meson P, containing quarks of flavors p and q, is given by 



If one could compute fp (i.e., fx or f D or f B ) with a reliably estimated uncertainty, 
a measurement of the partial width is tantamount to a measurement of \V pq \. 



B, D, K & CKM 
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TABLE 1. How to combine exclusive experimental measurements with calculations in 
(lattice) QCD to obtain the Cabibbo-Kobayashi-Maskawa matrix. 
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Similarly, the differential decay rate of a semileptonic decay is given by 
dTp^Hiv 



dq 2 



known 
factors 
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(2) 



where q is momentum carried off by the leptons, and T is the appropriate form 
factor. These processes do not suffer the helicity suppression of the leptonic decays, 
so the statistical error of the experimental measurements is smaller. Furthermore, 
lattice QCD can provide the q 2 dependence, at least when the momentum of the 
daughter hadron is not too large. If theory and experiment exhibit the same shape 
as a function of q 2 , one's confidence in the systematics increases qualitatively. 

Neutral meson mixing reveals a glimpse of the third row of the CKM matrix. 
For example, the mass difference of neutral B mesons is given by 
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B a fB a B B q \V t *V tb \ 2 , 



(3) 



m. 



where the flavor q can be either down or strange. The notation employing 
f B , and the "bag parameter" Bb q is historical and is taken from the kaon system. 
(In the so-called vacuum saturation approximation, Bb = 1.) Nevertheless, this 
formula, more so than leptonic B decay, motivates the interest in fs d and fs B - 

Figure 1 shows a time- line of calculations of /# with lattice QCD [10-30]. The 
results included in Fig. 1 have been selected with two criteria: Conference pro- 
ceedings, which are almost always followed (eventually) by papers in refereed jour- 
nals, have been omitted. Otherwise, I have taken papers whose authors were self- 
confident enough to quote a result in the abstract. The second criterion is not 

4 



500 



400 



300 



m 

^ 200 




100 



10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 

Ref. 

FIGURE 1. Time-line of calculations of fs with lattice QCD. Methods shown are extrapolation 
from m,Q < m c to nib (inverted triangles), the static limit uiq — > oo (squares), interpolation 
between m c and oo (circles), NRQCD (triangles), and that of Ref. [6]. Entries with a thick 
tick-mark on the horizontal axis control lattice spacing effects, as explained in the text. 



necessarily fair to cautious innovators, but, on the other hand, it generates the 
picture that is seen by outsiders. Figure 1 does contain a few exceptions to the 
second criterion, to include calculations that offered new technical developments. 

One can divide the time-line into infancy [10-14], childhood [15-22], youth [23- 
26], and adulthood [27-30]. (A similar, more discriminating classification has been 
made by Bernard [31].) In adulthood, with refereed publications starting in Novem- 
ber 1997, the scatter that characterizes the field through its youth has settled down. 
The results with a thick tick-mark on the horizontal axis use several lattice spac- 
ings [16,22,27,28,30] or a fully consistent implementation of lattice NRQCD [29]. 
Thus, one could say that Refs. [16,22] are mature results of the static approximation 
but, because the contribution of order l/m b is not negligible, not of 



HEAVY QUARKS AND LATTICE FIELD THEORY 

In the introduction, I explained that physical I/tuq effects and artificial a effects 
are intertwined and that a better appreciation of the intertwining was required 
before calculations of could mature. In particular, a theoretical analysis [6] of 
Wilson quarks away from the small-mass limit was needed to obtain the results in 
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Refs. [27,28,30]. This section summarizes some of the main ideas by comparing and 
contrasting the effective theories in the continuum and on the lattice. 
The effective Hamiltonian of QCD for heavy quarks can be written 

H = Mi** + *7 A * + */i*, (4) 

where * is the fermion field and h is given by an expansion in 1/mg, namely 

D 2 ill B {j-D,a-E} (D 2 ) 2 

~ ~2M 2 ~ ZB ^M^ + ^ SMI Zi ^Mf + " ' ' (5) 

corresponding to the kinetic energy, hyperfine splitting, spin-orbit splitting, rela- 
tivistic corrections, etc. 2 This result follows from a series of Foldy-Wouthuysen-Tani 
transformations, or from noticing the heavy-quark symmetries of QCD, writing 
down allowed operators with arbitrary coefficients, and matching physical observ- 
ables to standard QCD. With radiative corrections, one finds z = 1 + 0(g 2 ). 

The rest mass of a quark is M\ and the kinetic energy of a quark is p 2 /2M 2 . It 
is convenient to call M 2 the kinetic mass, even though Lorentz invariance implies 
M 2 = Mi. Let us write H = Mi** + Hq. The rest-mass term and Hq commute, 
even in the interacting theory, so eigenstates of H are simultaneously eigenstates 
of Hq. Thus, one can drop the rest-mass term or readjust Mi according to con- 
venience, without changing the dynamics of heavy-quark QCD. The physically 
relevant parameter is the kinetic mass, which one adjusts so that M 2 = tuq. 

The lattice Hamiltonian takes the same structure as in Eqs. (4) and (5), but 
with changes to the operators' coefficients. One can express this by replacing h in 
Eq. (4) with h\ at = h + Sh and writing 

5h = aftfliE ■ B + a\. D .{7 D,a E} + a%(D 2 ) 2 + a 3 w A £ D\ + • ■ ■ , (6) 

i 

where the coefficients b = b(niQa,g 2 ) and w = w(rriQa, g 2 ) are functions of the 
(lattice) quark mass m^o and the gauge coupling g 2 . The same operators appear, 
as well as others, such as the last one, that break rotational symmetry. Part of the 
craft of the numerical work is to adjust the underlying lattice action so that these 
artifact coefficients b and w vanish, at least to some accuracy. 

In lattice NRQCD and HQET this pattern arises by construction, in particular 
one sets Mi = 0. It is fairly straightforward to adjust the bs and ws to vanish at 
the tree level of perturbation theory in g 2 . Beyond the tree level it is still possible, 
but arduous perturbative calculations are needed. Power-law divergences appear 
in loop diagrams, so one cannot take a — > by brute force [2]. That means that 
the lattice artifacts can be removed only by further refinements of the NRQCD 
action. For heavy-light mesons, the action of the light quarks must be improved to 
a consistent level [32], as in Ref. [29]. 



2 ) Analogous expansions are introduced for the operators mediating electroweak transitions. 
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With actions for Wilson quarks, such as the clover action 3 or the Wilson ac- 
tion itself, the pattern sketched here is less immediate. It is guaranteed, however, 
by the heavy-quark symmetry of the lattice action. Because the lattice violates 
Lorentz invariance, the rest and kinetic masses differ; in practice Mi < M2. The bs 
and ws are, however, bounded for all masses. As rrtQa — > 0, these coefficient func- 
tions go to a constant 4 or vanish like (rriQa) S0 , for some integer sq, by Symanzik's 
standard analysis of cutoff effects. On the other hand, as rriQa — > 00, they vanish 
like (l/mga) s °°, for some Sqo, by heavy-quark symmetry [6]. The full functional 
dependence on mqa is calculable in perturbation theory and, sometimes, nonper- 
turbatively. The coefficient functions bs and ws depend on the details of the lattice 
action. For example, with the Wilson action bs 7^ 0, but with the clover action one 
can adjust an unphysical coupling until bs vanishes. 

LEPTONIC DECAYS 

Let us now return to Fig. 1. The plotting symbols distinguish methods. Squares 
denote calculations in the static approximation, and triangles denote calculations 
with lattice NRQCD. Inverted triangles [10,11], circles, and diamonds denote treat- 
ments of numerical data from the Wilson or clover action. 

Because the latter connects naturally with the effective theories, Ref. [6] sug- 
gested treating the bottom quark on the lattice by adjusting the bare mass until 
M2 = mft. In addition, a suitably normalized operator, essentially one built from 
the heavy-quark field must be used for the current. Results with this approach 
are denoted with diamonds. The discretization errors are of order Aqcd« to some 
power, from matrix elements of the operators in Eq. (6). (These effects are then 
multiplied by a coefficient b or w, which is a number of order 1 or, in limiting cases, 
smaller still.) With this interpretation of the numerical data from the clover action, 
the lattice results for f B are nearly independent of the lattice spacing [27,28]. 

In earlier work with Wilson quarks, the tuning of the mass and the normaliza- 
tion of the current introduced discretization effects of order rriQa or (rriQa) 2 . The 
authors minimized these lattice artifacts by reducing the quark mass of the lattice 
calculations below that of the charmed quark. This is still large: m c a ~ 5Aqcd«- 
Then the results were extrapolated up to m^, with fits to l/mq expansions, either 
with (circles) or without (inverted triangles) the help of the static value. This in- 
tertwines lattice artifacts and 1/rriQ effects in ways that depend on details of the 
fits, not least because one cannot verify whether the 1/rriQ expansion converges for 
the low quark masses, on which the fits are based. 

Table 2 tabulates the mature results [27-30] for f B , f Bs , f D , and fo s , along with 
an average [34] of results from experimental measurements of T^^^u, combined 
with \Vcs\- While it is tempting to average the results, it is subtle to do so properly, 



3 ' The clover action adds a term to the Wilson action to eliminate the leading lattice artifact [33] . 
4 ) If the limit is a constant, the lattice artifact still vanishes, owing to the explicit powers of a. 
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TABLE 2. Compendium of recent results for decay constants of heavy- light 
mesons. The experimental average comes, in fact, from measurements of |V^ s |/£>, , 
which then take \V CS \ from unitarity or from neutrino production of charm 
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because the systematic errors are largely, but not entirely, common. If you must 
have an average, consult Draper's review [9], or use your own eye. 

The central values in Table 2 (except the experiment!) are in the quenched ap- 
proximation. The error bars, however, reflect estimates of the associated uncer- 
tainty. The best estimation is that of the MILC Collaboration [30] who have some 
lattice results including up and down quark loops. (The strange and heavier quarks 
are still quenched.) The partially unquenched data sometimes lie higher than the 
quenched data, sometimes not, leading to the very asymmetric error estimate. 
These results are encouraging, not least because they suggest that the wait for a 
fully unquenched calculation will not be too much longer. 



NEUTRAL MESON MIXING 

Because the calculation is technically more demanding, the literature contains 
fewer calculations of the bag parameters of B% and B® meson mixing than of the 
decay constants. Recent publications report Bg s /BB d ~ 1 and 

1.03 ±0.06 ±0.18 [35] 

1.40 ±0.06 [36] ( ] 

1.23 ±0.05 ±0.15 [36,35] ' y J 

1.17 ±0.09 ±0.05 [37] 

where Bb a is the renormalization-scheme independent combination. The first three 
entries use the static approximation. The third entry comes from an analysis in 
Ref. [36] of the data in Ref. [35]. The last entry uses lattice NRQCD and finds, 
additionally, that the dependence on I/tuq is not large. See Ref. [9] for more 
results. There are also results for the ratio of matrix elements: 
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\m 2 B J 2 B B Ba f 1.38 ± 0.07 [35] 

lmlJl d B Bd \ 1.76 ±0.10 IS [38] ' 1 } 

On the whole, the impression is that more work needs to be done to gain control 
over the systematic errors; calculations of B B are not as mature as those of fg. 

The kaon's bag parameter Bk is needed to predict a measure of indirect CP 
violation in K — > tttt. Given Bk from (lattice) QCD, a measurement of Ek traces 
a hyperbola in the complex Vta plane. The history of Bk has similarities to that 
of Jb'- numerical work was a greater challenge than initially hoped, and theoretical 
insight is needed as a guide. For example, the lattice-spacing dependence (with 
staggered fermions) is surprisingly steep, but one now knows to extrapolate to the 
continuum limit in a 2 (rather than a) [39]. Two recent calculations find 

B fNDR 2CeVl-J°- 62 ±0.02 ±0.02 [40] 

iMNDR, 2 GeV) - j 628±0 042 [41] , (9) 

where the main uncertainty comes from the continuum extrapolation. Quenching 
and degenerate quark-mass effects may each lead to underestimates of order 5% [42]. 

SEMILEPTONIC DECAYS 

Semileptonic decays are wonderful for learning about the first two rows of the 
CKM matrix. Because a hadron is in the final state, lattice calculations of the form 
factors are more difficult than decay constants, but not much more difficult. 

In these decays there is an additional kinematic variable, the momentum \p'\ 
of the daughter hadron (in the parent's rest frame). The Lorentz invariant q 2 is 
linearly related to \p'\. Until recently, calculations of these form factors were done 
with mQ^m c and extrapolated up to m c or with 1/mg expansions. Since 
the kinematically allowed range of q 2 depends on the heavy-quark mass, another 
extrapolation is made. It is clear that the extrapolations once again intertwine 
artificial a effects with physical l/mg and q 2 dependence. The details of the inter- 
twining are not transparent (to me, anyway), so I reserve comment and direct the 
reader to reviews by Onogi [7] and by Draper [9]. 

More recently, calculations have been done with lattice NRQCD or with Wilson 
quarks interpreted [6] as suggested by Eqs. (4)-(6). These techniques are especially 
powerful here, because the ability to compute directly at the B mass decouples 
extrapolations in q 2 from the mass dependence. Calculations are underway for 
form factors of heavy-to-light transitions, such as D — > irlv and D — > Klu, which 
together yield |F CC ;/F CS | [43,44], and B -> irlu, which yields \V ub \ [43,45,46]. 

I would like to conclude with a preliminary result [47] on the zero-recoil form 
factor for the decay B — > Dlv, which will improve the determination of |V^|. A 
similar study of B — > D*lu is in progress. Until now there have been calculations 
of the shape of the form factors for B — > D^*Hu, but the normalization has been 

9 



computed only poorly, in perturbation theory It is possible, however, to handle 
almost all of the normalization nonperturbatively. For example [48], 



(D\V \B)(B\V \D) = h B + - D {l)h D + - B {\) = ^ 
(D\V \D)(B\Vo\B) h^ D (l)h^ B (l) 1 + 



/i^(l)| 2 , (10) 



and here at Fermilab we have found analogous ratios for h B ~^ D (l) and h^~* D *(T). 
For the class of actions considered in Ref. [6], the remaining radiative corrections 
have been computed [49], and they are small. 

In Eq. (2) for B — > Dlv one requires T = h + — (tub — ra£>)h_ / '(mg + mo). 
The advantage of the new method is that, in effect, it calculates not ^"(1) but the 
deviation of .F(l) from 1. We find 

= 1.069 ± 0.008 ± 0.002 ± 0.025, (11) 

where the uncertainties are Monte Carlo statistics, tuning of the quark masses, and 
a certain parametrically small contribution omitted from /i_. Uncertainties from 
lattice artifacts and from quenching have not yet been taken into account. 
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